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Thermal Diffuse x-Ray Scattering in Simple Metals * 

by 

David M. Straus and N. W. Ashcroft 
Laboratory of Atomic and Solid State Physics 
and Materials Science Center 
Cornell University, Ithaca, N.Y. 14853 


Abstract 

Calculations are reported for the ionic structure factor and x-ray scat- 
tering cross section of sodium (at T=0°K and 90°K) and lithium (both isotopes 
it T=0°K) within the harmonic approximation. An evaluation of the appropriate 
lisplacement- displacement correlation function by the special point method 
lircumvents the need for a multiphonon expansion. In the case of sodium, the 
structure in the one-phonon scattering is straightforwardly accounted for, and 
m approximate expansion is obtained for all multiphonon scattering. By treat- 
ing core and conduction electrons on an equal footing it is shown that informa- 
tion on the conduction electron system is present in the forward scattering 
component . In lithium the one-phonon cross section at small angles aids in 
the determination of the effective electron-ion interaction. 
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I INTRODUCTION 

For some years x-ray thermal diffuse scattering (TDS) has been used as a 

1-4 

probe of lattice dynamics in simple materials. • Although information on the 

phonon frequencies and polarizations (and also the extent of anharmonicity) is 

4 5 6 

contained in the TDS, ’ it is generally hard to extract. The cross section 

for the scattering of x rays intimately involves the static structure factor 

of the ions, S ion GO. The purpose of this paper is to present calculations 

of (i) ST on (k); a nd (ii) the x-ray scattering cross section for Na and Li in 

the harmonic approximation and in their ground states. The significant features 

8 9 

of the calculation are the use of a special point technique * in the compu- 
tation of the equal time displacement-displacement correlation function <^u ^ 
(which enters into S ioir (k) ) and the separation of the scattering cross sec- 
tion into contributions from core and valence electrons. In particular, the 

4 

special point technique enables us to avoid the customary expansion of the 
inelastic part of S ion- GO into terms involving the 'scattering of a definite 
number of phonons. We determine the "one-phonon" term explicitly, but we can 
also calculate all higher order processes without recourse to expansion. Fur- 
ther, our treatment of the contribution of the valence electrons to the cross 
section shows that x-ray scattering should yield information, in light metals, 
on the effective electron-ion interaction, as we demonstrate for the partic- 
ular case of Li. 

Section II contains a derivation of the x-ray scattering cross section 
da/dQ in a model of a simple metal which distinguishes between bound and con- 
duction electrons. In Sec. Ill we outline the calculation of S ^ 0 n Gj^ using 
the special point technique (discussed in detail in the Appendix) , and compare 
it with the other non-expansion techniques in the literature. Section IV pre- 
sents numerical results for S ion (k) and do/dG for Na (at two temperatures) 
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and for both isotopes of Li . We draw particular attention to the secondary 
maxima associated with the one-phonon term as observed in certain crystallo- 
graphic directions. These maxima have special importance in the determination 
of the electron-ion interaction, of Li, and also give information about specific 
portions of the phonon spectrum directly. 
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II. THEORY 


The differential cross section for scattering of a photon from a solid of N 
ions in volume V (at T = 0°K) is proportional to the space-time Fourier transform 
of the Van Hove correlation function G (r.t): 

4 0 v ~ ' 

2 

d o C f, 3 P 00 

dfidm " V J d r J dt G e <S>t)exp[ik-r - imt] , (2.1) 

— CO 

where C is a constant,^ ^ 


G e (r } t) 


d^x <(ft(x,0)h.(x + r,t)^>. 


( 2 . 2 ) 


and 


~ = ~i “ ~f ; w = io i - u f . (2.3) 

We are considering the cross section per unit volume for scattering a photon 
of momentum Tik^ and energy hoK into a solid angle dW with energy loss between 
“hw and h(w + dm) . The quantities hk^ and 'hoo^ are respectively the momentum and 
energy of the scattered photon. In Eq. (2.2) n(r,t) is the total electron number 
density operator and the brackets^ prefer to a ground state average. Introducing 
spatial Fourier transforms, 


d 2 o 

dfidu) 


n c 

_ c. 

v -0 


I 

dt e _lait <n(-k,0)n(k,t)^> , 


(2.4) 


-where n(k) is the Fourier transform of n(r) . 

We separate n(r) into contributions from core and valence electrons, and we 
treat the core electrons as if they were rigidly attached to the ions. Any core 
excitations or distortions of the ions are therefore neglected; should these occur 
they must be calculated separately. In practical terms this means that in comparing 
experiment and theory the Compton scattering from the core electrons must first be 
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subtracted from the data. In addition we invoke the adiabatic approximation, so thal 
the conduction electrons (ce) are always in a ground state appropriate to an in- 
stantaneous' ion configuration (ion) . By virtue of the rigid ion approximation we 
may write: 


h(k,t) = Fe ife ‘?i (t) f(k) +ft (k,t). (2.5) 

x ** ce ~ 

Here f(k) is the Fourier transform of the average core electron density about a 
nucleus at the origin, and R. (t) refers to the instantaneous position of the ion 
labelled i. From Eqs. (2.4), (2.5) and the adiabatic approximation, we then find: 


d 2 a 
d&dw 


% - I “dte-^C ( 

^ ^ ion 




+ 


+ 


<(e f(-k) <^h 

i ' 

J_„ dt 


(k,t)\ V 

ce ~ / ce /ion 

(-k,0)\ V ) 

- / ce/ ion 

n (k,t)\ 
ce ~ / ce,ion 


( 2 . 6 ) 


We suppose that the interaction between conduction electrons and ions can be 
represented by a weak pseudopotential with Fourier transform v(k) , (as is the case 
for many simple metals) . The density response may then be calculated to linear 
order in v(k): 


with 



(k,t) 



X 1 (k) v Qs> 


Z e i ^?i (t) , 


(2.7) 


X x (k) 



( 2 . 8 ) 
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G(k,0) being the static dielectric function of the uniform interacting electron 
13 

gas. Equations (2.6), (2.7) and (2.8) now give: 


d 2 o < V 
dSldw C 


J 


-iwt 


R- i “‘ ( ~ l( ° ) ( l f (k)| 2 + 2£(k) Xl (k)v(k)) 

L- ±j 

+ (ft (-k,0)ft (k, t) ) . 1 

\ ce' ce - V ce, ion]. 


(2.9) 


In a typical x-ray experiment all the radiation emerging at a given angle 
. ... 14 

is initially measured. All possible energy transfers (on the scale of typical 

electron and phonon energies) hoJ are therefore included, and we pass from the cross- 

2 

section for energy loss fiw (d a/dfidw) to the total angular cross section (do/dfi) : 

da - O" fair ■ ft 1 e " 1 -' ion < |f (k) I 2 + 2f (k)x 1 (k)v(k) ) 


+ ( ft (-k) n (k)^ 

\ ce J ce ~ / ce,i 


ion ] 


( 2 . 10 ) 


Note that the last term is usually considered part of the Compton scattering, and 

16 

is therefore generally subtracted from the primary data. What will become appareni 
in Sec. IV, is that the value of the last term in Eq. (2.10) (the valence electron 
correlation function) should be readily obtainable from x-ray measurements. The 
theoretical results we present are therefore best compared to data from which 
only the ionic Compton scattering has been subtracted. 

The last term in Eq. (2.10) is difficult to calculate for interacting electrons 

in the presence of the ions. For purposes of illustration we use the free-elecfron 

, 13,17 

value 


N S (k) 
e e _ 


* ( fi ce < -W fl c e W>ce. 


free , 


s e 0s) 


3k 


1 k~ 


4k f 16 k f " 3 1 ’ 


0 < k £ 2k f ; 


( 2 . 11 ) 


= 1, 


k> 2k f< 
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Here N e is the number of electrons, and kf the Fermi wave vector. Setting 
(for-a monovalent system) the number of electrons N e equal to the number of 
ions N, Eq.s (2.10) and (2.11) give, us the final result: 

* 55 • S *2Sc “ S ion<« • <l f< « I 2 + + S e<S>’ < 2 ‘ 12 > 

where we have set 

S. (k) e _1 ^* . for k 4 0, ' (2.13) 

ion ~ N V j ion ~ 

It should be clear that except for the elements of lowest atomic number 
(e.g. Li), S„(k) makes a small contribution to W for all but the smallest wave 
vectors k < 2k f . 



Ill CALCULATION OF IONIC STRUCTURE FACTOR 


We now proceed to a calculation of S. (k) in a model in which the solid 

ion ~ 

is treated as a harmonic crystal. Letting R i = + u^, where is the 

equilibrium position of the i* k ion and u; its displacement, 


<'-i 


S. (k) = i £ e-^Sij ~ ^j>> 

ion ~ N ' ' ion 


(3.1) 


Here, , = X. - X., and the average is to be taken over the states appro- 


priate to a harmonic crystal. With the definitions 
< s5 u i U j^a ^ U i “ U j' ) ^ ion ~ ( ~i ~ 


and 


we 


S." <U 1’ V V- 

4 19 

have the result ’ : 


- SJ>) = e-K'fWi-Jl' 

\ f ion * 


where 


ap ~i 


h 1 

“ ?. 4 ) = jjjj £ (l-cosq*X_)e /v (qj)e a (qp) v coth [|{3hju (qj ) ] ,. 


BZ 

£ 

qj 


(3.2) 


(3.3) 


~ ^ij a 


U) (qj) 


P = 1/^T, 


(3.4) 


and M is the mass of an ion. In Eq.(3.4) 0) (qj) and e(qj) are the frequency 
and polarization vector of the normal mode of wave vector q and polarization 
index j, (j=l,2,3). The q sum extends over the entire first Brillouin zone. 
Using the translational symmetry of the lattice, Eqs. (3.1) -(3. 4) yield: 

S (k) » £ e -i ~ Si e _ ^ k P Xa 3 ( 5. i:) . (3,5) 

ion ~ 


8 
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Next we separate as follows: 


WS’-W 05 - v^’ 


with 


(3.6) 


A ae ( 5i 5 = -w q E . Vs^ysh ;^|jy ooth gsimsih <=°= <a-&> 


Note that 


V ( 2iJ ) 2 ^ia U je) ion. 


(3.7) 


(3.8) 


We see therefore that A Q (X) is the displacement-displacement correlation func- 

CXp ~ 

tion. for two ions separated (on average) by X. Clearly A „(0) is the dis- 

~ ap 

placement-displacement autocorrelation function. For a cubic system, 

ion = 6 af? 5 4 'i) • <3 ' 9) 


so that 


W 0) 


1 ii 


§ a£ 3 MN q j u) (q j ) 


coth[^phu (qj) ] = 6„ c A° 

(Xp 


(3,10) 


o 4 *~2W 

This defines A , which is closely related to the Debye-Waller factor e : 


av = VeW 0) = 


(3.11) 


Substituting Eqs. (3.6)-(3.11) into Eq.(3.5), we find: 

(k)= !,■%'!» 


ion- ~ 


= S e~^ 2h '° e +ik a k p A aP ( 2l ) 


(3.12) 


To proceed from this point the usual approach is to expand the last 

exponential in a power series in A _(X.). The leading (i.e. constant) term 

<Xp ~ a - 
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gives the elastic (Bragg) scattering peaks, the second gives the one-phonon 

scattering, the third the two-phonon scattering, and so forth. Beyond the one- 

phonon contribution each term is increasingly laborious to evaluate. We can 

avoid this expansion, however, by writing S. (k) as follows: 

ion 


S. (k) = S e ~ ~i e 

ion w i 


-ik-x, -£k 2 A° + z g-ik-X* e ik A° [|k^k o A^ fi (Xi)] 

i 


a p ap 


+ s e - 4k2A ° -1 -i* lUL a <Jl>] 


= S (k) + S (k) + S (k). 

O 1 r*t M rv 


'a 

(3.13) 


Here S (k) gives the elastic scattering, 

l k 2 A° 

i.e. S (k) = Ne 2 £ v , 

o k f *“■ 

r*-' rs^ 


(3.14) 


the K being the vectors of the reciprocal lattice. The one-phonon scattering 

term S_ (>k) is easily seen to be: 

* X 

S (k) = e" 2k A ° — k k 2 e [q(k)j] e [q(k)j]~ Xpppj coth[|hpa)[q(k) j ]1, 
i v x/ 2M a a ( 3 <^~ u)[q(k)j j l - 1 

(3.15) 

where q (k) is the vector k reduced by an appropriate K to the first Brillouin 
zone (i.e., q(k) = k - K) . Finally the remainder S (k) will be calculated by 
direct computation of AlB(Xi), so that all higher order phonon terms are auto- 

- /v 

matically taken into account. The reason for adopting this procedure is to 

/ 

assure convergence in the sum over i in S (k) . This will be clarified in what 
follows . 

Our method of calculation of A „(FU) and A makes use of special points 

ap ~ L 

8 9 

in the first Brillouin zone ’ to evaluate the integral of Eq.(3.7). By 
calculating the integrand at these relatively few special points, one obtains a 


good approximation to the entire integral. This procedure differs markedly 
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from ordinary numerical integration in that (as shown in the Appendix) one 
is effectively using an expansion of .the integrand in symmetrized plane waves. 


In connection with this method we draw attention to the behavior of V^^(X) 
for large X.- At large X the dominant contribution to the integral in Eq.(3,7) 

r -' /v 

20 

comes from small q, and it can be shown that at T=0 °k, 


Lim 

X-*» 




(3.16) 


Thus to ensure convergence in S,„(k) it is necessary to make the separation 

M <■- 

.indicated in Eq.(3.13). 

The method may be compared with the non-expansion calculations of S^ Qn (k) 

21 

by (1) Lomer, who calculates the ionic structure factor directly using the 

22 

results of a computer experiment; (2) Semenovskaya and Umanskii, who cal- 
culate A „(X) in closed form for a model sinusoidal phonon dispersion law; and 
Ctp 

23 

(3) Reid and 1 Smith, who calculate the multi-phonon scattering S„(k) for 

M r** 

crystal's whose sizes range between 100 and 1000 unit cells. Their- evaluation 

of A (X) is acheived by summing over only those q corresponding to the normal 

ap ~ ~ 

modes of such a finite crystal. By separately calculating the (£~*0 portion 

of the integral in Eq.(3.7)., they find that a crystal of 500 unit cells gives 

essentially the same S (k) as an infinite crystal, for q(k) belonging to the 

set of normal modes of the finite crystal. 

The method of Reid and Smith appears to be the most accurate and practical, 

but has the disadvantages that one can calculate S^(k) at relatively few points, 

and that the matrices A „Q0 for a real crystal are inaccessible. We are able 

Ctp ~ 

to circumvent these limitations by directly calculating the correlation matrices 

A _ (X) - (These are of considerable interest, of course, in a wide range of prob- 
ap ~ 5 

lems . ) 

We illustrate, the method by its application to Na and Li. In both cases the 
phonon spectrum was calculated from a force-constant model designed to fit the 
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experimental data. The corresponding S. (k) has "been calculated for Na at 

ion ~ 

two temperatures (0°K, 90°K) and for both isotopes of Li (at T—0°K) . 

In the case of Na the force constants were those that fit the data at 

24 25 

T=90°K . A simple estimate (supported by some theoretical results ) in- 

dicates that the change in phonon frequencies between 0°K and 90°K is everywhere 
less than the experimental error. Hence the only effect of temperature we 


allow is through the hyperbolic cotangent function in Eq.(3.7). 


To simplify 


the calculation we use the T=0°K value of A .(X.) for X.^0 in the 90°K cal- 

o£ ~1 

27 

culation, but use the t= 90°K value of A (0). - The 90°K 

CC.p - 

resultsare therefore meant to be indicative of the effects of temperature, but 


they are only approximate . We use the value of r g determined from the 5 °k 

28 

lattice constant measurement i.e. r = 3.931 a.u . ( r is defined by 

s s 

-rr(r a ) 3 = — , where a n is the Bohr radius.) 

s °N’ ° 

fci 

7 

The force constants for Li were similarly taken to be those which fit the 

29 

experimental phonon dispersion measured at T = 98°K.The values of r was also deduct 

s - 

28 

from the lattice constant, in this case at 78°K (r = 3.248 a-.u.). To cal- 

' s 

culate S. (k) we have set T=0°K, In order to obtain S, (k) , A°, and A ^(X.) 

ion ~ 1 ~ ap ~i 

6 

for Li, we have assumed that both substances are truly harmonic. This gives: 

-I 

(ycc M 2 , 

X 

A n (X) o: M 2 for all X 
ctg ~ ~ > 


and S. (k) <= M 2 . 
1 ^ 


(3.17) 




IV RESULTS 


In this section we present numerical results for both S (lc) and the 

ron ~ 

x-ray scattering cross sections for Na and Li. The structure factor calcula- 
tions were carried out as described above. As regards the cross sections, we 
give two sets of results. One corresponds to the theory outlined 1 , in Sec. IP: 
da V 1 


W = 


dQ N 


2rrC 


SCk)( jf(k)| 2 + 2f(k)v(k) v(k)) + S (k) , 


(4.1) 


while the other corresponds to the more commonly used expression: 

W e te J -i-) = S(k) |f (k) f 2 

a \dQ N 2 ttC/ 0 ~ 1 a ~ 


(4.2) 


Here f (k) is the Fourier transform of the average electron density of an 
a ~ 

assumed neutral atom, and we write (and shall continue to do so) S(k) in place 

of S. (k) . Both the ionic (Na , Li ) and the atomic (Na, Li) form factors 
ion * — * 

31 

were taken from Ref. 30. The Geldart and Vosko modified form of the Hubbard 
dielectric function e(k) was used, as well as an empty core pseudopotential 

32 

to represent the effective electron-ion interaction. Figures 1-7 show 
S(k) for Na, and Figs, 8-10 show S(k) for both isotopes of Li. Numerical 
results are listed in Tables I-IV. We present both cross sections V and W 

a 

for Na (at T = 0°K) in ' Figs. . 12-15, and in Figs. 16-18 we show V for Li 
(at T = 0°K) with two choices of the core radius appearing in the empty 
core pseudopotential. 

The most noticeable feature of the structure factor plots is the sizeable 

structure between the Bragg peaks along all directions except the [100] and 

[110] directions (for a bcc lattice) . These maxima are a direct consequence 

33 

of the behavior of the one-phonon term. Their occurrence is completely 

34 

general, and has been noted for quite some time. For the sake of simplicity, 
however, we can most easily explain them in terms of a (polarization-independent) 


13 
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Debye model. Here (at T = 0°K), 


S..(10 

X 


_4k 2 A° _h k 2 i 

2M W(a(k)) ’ 


(4.3) 


where uuGq,) = cq is independent of polarization and c is the approximate speed 

of sound. We have plotted in Fig. 11 lines along which the function 

l/ c lfl(h)| has constant value for a (001) plane of the reciprocal lattice of 

a bcc crystal. In any direction (except [100] .and [110]), and as a consequence 

of periodicity alone, the one-phonon term displays secondary maxima as one 

passes over the ridges of the function shown. Replacing Eq. (4.3) with 

Eq. (3.15) introduces three frequencies (one for each polarization j) at 

2 

every point, each weighted by the factor [k* < e( j&(k) ) ] . For example, Fig. 11 

37 

would indicate two secondary maxima between the points (1,0,0) and (3,1,0), 

whereas Fig. 7 shows only one. The value of the one-phonon term at the point 

along [310] marked P on Figure 11 is determined by the phonons at the point 

& = (i,^, 0) in the first Brillouin zone. At (^,4,0) Na has an anomalously 

24 

low transverse frequency. Furthermore, since q, is nearly perpendicular to 

2 

the [310] direction, the factor £ [e(jq,)*k] will select out the transverse 

3 ~ ~ 

frequencies. The resulting single, large, maximum swamps any other effects. 

Thus we see that any particularly low phonon frequency will cause a sequence 

of one-phonon maxima along the appropriate direction . This property of the 

• 33 

one-phonon scattering has been widely used to study soft modes, but the 
discussion is often set in real space. In terms of identifying the maximum 
with a particular vibrational mode we see that it is advantageous to treat 
the problem in reciprocal space. 

The comparison of W and W for Na in Figs. 12-14 shows that at large 

s 

k the only significant difference is a shift arising from the term S (k) 

^ e ~ 

in W, which is a constant for k > 2k f . However, at small £ Fig. 15 shows 
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that the presence of S (k) in W contributes to a difference in shape between 

0 

W and W . The small k portion of the x-ray cross section (with only ionic 

3. 

Compton scattering subtracted out) thus gives us information about the conduction 
39 

electrons. Note also that for Na the presence of the pseudopotential v(k) in W 

seems to make little difference in the final cross section. This is not so 

for elements of very low atomic number. For example, in Figs. 16 and 17 
7 

we plot W for Li at low values of k, for two choices of the core radius ap~ 

40 

•pearing in the empty-core pseudopotential. The maximum percentage difference 

is slight...in;.both cases, but in Fig. 17 the actual shape of the one-phonon 

maximum is noticeably altered. In fact, the differences between pseudopotentials 

will always be most noticeable in low k one-phonon maxima. In order for v(k) 

to have any influence in Eq. (4.1), we need to 'have k < 2k_ (otherwise y_ is 

f 1 

exceedingly small) and S(k) to be not too small. Figure 18 emphasizes this 
point: Here we plot W - S e (k) , so we subtract all the Compton scattering. 

What remains shows a marked dependence on the pseudopotential. 

We should discuss the relative composition of the TDS (i.e. of S(k)). 

/v 

Figure 1 shows the contribution of the one-phonon term, and we see that at 
large k the many-phonon terms become quite important. From Eqs. (3.11) and 

ro 

41 

(3.13) : 

S M (k) = e "^ A tXp (0) [e^ k PV (0) - 1 - £k a k p A af3 (0) 


+ vg Q e e 2 ka k (3^(Xf3 e 2 ^a k pAxg - 1 - Jk^kgA^ Q^) J . (4.4) 


From the Appendix we also note that for Na, Tr/^Qj;^) < < TrA^ (0) . (f or ^ 4 0>. 
Typically at least 90% of S,,(k) in Na comes from the first term .in Eq. (4.4), 


S M (k) 

M ~ 


PS 


1 - e 


i.e . ; 




(4.5) . 
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In Eq. (4.5) we have confirmed a well known approximation. 

In spite of the fact that the X- sum in s M (k) converges roughly as 
4 

E. 1/jr . , we have found it adequate to take only 9 shells (136 vectors) 
in the sum. (Taking only 7 shells changes S(k) for Na by considerably less 
than 1%, for example.) This can be understood by noting that 



(4.6) 


where X and X are typical vectors in the first and ninth shells. The point 

2 

is- that the asymptotic limit of A (X,. ) (“ 1/X . ) is only reached at large 

CCp i 

where the structure factor is almost independent of the contribution of 

the remaining shells. In addition, the sum actually converges more quickly 
4 -ik*'Xj 

than E-; 1/X. since the term e in Eq. (4.4) introduces (except for k = K) 


.considerable self -cancellation. 



V DISCUSSION 


The extension of our method of calculation of the ionic structure factor 

to systems without cubic symmetry and to systems with a basis is completely 

9 

straight! orward . (Special points have been found for systems of hexagonal 

symmetry, and they can be generated for systems of any symmetry.) The 

occurrence of one-phonon maxima is equally general. The ability to calculate 

the A _(X.) by a procedure which avoids a difficult three-dimensional numerical 
Ctp '■'tL 

integration should prove valuable in a variety of contexts, including, for 

19 

example, the self-consistent harmonic theory of phonons and the computation 

42 

of static lattice Green's functions. 

Much of the theory of x-ray scattering from simple metals presented in 

44 

Section II can be extended to liquid metals. Egelstaff, March and McGill 
have derived a formula for the x-ray cross section in liquid metals that is 
identical -to -Eq-. -(.2.6), except that they do not make the adiabatic approximation 
in the terms involving the correlation of conduction electrons with the ions. 
Making that approximation, and introducing the pseudopotential v(k) , we con- 
clude that Eq. (2.12) is as valid for liquid metals as it is for crystals. 

Finally, our calculation has neglected possible anharmonic effects. Those 

19 

anharmonic terms which are retained in the self-consistent phonon theory are 
in a sense taken into account hez’e. The formalism we have presented is not altered 
by zising the self-consistent theory, but the frequencies are changed from 

OK 

their harmonic values. In the case of sodium, this change is small. Other 
anharmonic effects are not taken into account. For example, the interference 
between one- and two-phonon scattering can cause a noticeable change^ 5 ^ in 

n 

s ^ on (k) . As shown by Glyde, however, it amounts to .only a small shift in 
the one-phonon scattering for Na at high temperatures. Since both the an- 
harmonic frequency shifts and the inverse phonon lifetimes become quite 


17 
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small at low temperatures, the size of this contribution should decrease 
correspondingly. Interference effects, as well as other effects due to an- 
harmonicity, may of course be of somewhat greater importance in the case of 
lithium. 



APPENDIX 


8 9 

We briefly review the special point method, ’ which was designed for 
the integration of quantities varying slowly over the first Brillouin zone (BZ) . 
Here, by a slight modification, we use it to evaluate the integral of oscillating 
•functions (see Eq. (3.7) ). 

The general integral to be evaluated is: 

1 bz r 8 

n 5 f< 9> ='(iff)vjBz “ * <A -» 


where f(a) is assumed to be invariant under the operations of the crystal 
point group, and Q c is the primitive cell volume. (If f(&) is not symmetric, 
it can, of course, be easily symmetrized.) One expands f(ft) in symmetrized 
plane waves A m (&) : 

OD 

f(&) = f + XL f A (q), (A. 2) 

™ o m=l m m ~ 

with 

A (a) = S e 1 ^’^ (A. 3) 

m ** ^;m 

and 


f 

m 


i r 

N m ^ BZ 


d 3 q f(&) A m (a). 


(A. 4) 


X;m refers- to all lattice vectors ^ with the same length that are related 

t 

by point group operations. is the number of vectors in this m shell, 
and the sum in Eq. (A. 2) is ordered so that those shells with lowest X come 
first. 

A set {a±} of special points is defined as a set of n points in the BZ 
with associated weights which satisfy: 
n 

.XL a. A (q.= ) = 0 for m = 1, . . . , N, (A. 5) 

1=1 l ni ~- L 
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n 


a. = l. 

1=1 i 


Using Eqs. (A. 5) and (A. 6) in Eq. (A. 2): 


n n 

f = .E a. f (q . ) - .E.. a.A. T (q . ) f. r .. + 
o i=l i i=l i N+l ~i N+l 


(A. 6) 


(A. 7) 


Since f is the desired integral, Eq. (A. 7) gives an approximation to the 

integral consisting of an evaluation of f(q) at a (small) set of points. 

The first neglected term can be shown to be + f ^ ^ . Not all coefficients 

f for m > N have been neglected, as Eq. (A. 5) is always satisfied for an 

infinite number of shells. The index of the first shell for which Eq. (A. 5) 

is not satisfied is N+l. With increasing number of points n in the set, both 

the number and the magnitude of the neglected terms become smaller. 

At T = 0°K , TrA _(0) E l/ti)(jq) is a smooth -function, and we may 
CCp qj 

apply the special point method. Although the expansion coefficients f^ 

decrease slowly with increasing m for large m, they are much smaller than 

TrA^ (0) itself. Thus we expect increasing the number of special points n 

to have a small effect on TrA (0) • From Table A-I we see the convergence 

dp 

is more rapid for T = 0°K than for T = 90°K. 

The calculation of A „(X.), ^ 0 is more troublesome, and we illus- 

(Xp 

trate by examining the Trace of this matrix. Symmetrizing the integral of 




TrA^Xi) « 


BZ 

£ 

& 


l 

w(jq) 


A i 


(A. 8) 


Applying the special point method to this integral means neglecting some of the 

coefficients T whose form is (we are at T= 0 °k): 
m 

BZ 

i cc v — - — 

m fi)(jq) 

• ro 


A i^m * 


(A. 9) 
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Now A^A^ is itself a sum of symmetrized plane waves: 


A. (q) A (q) = £ a . (i > m)A (q) , 

~ m ~ -i j j 


(A. 10) 


where the first j for which a.(i,m)^0 is that for which X. = X - X. j. 

J ■ — j ^-m ~i 

From Eq.s (A.8)-(A. 10) it is clear that the T. for large m will be much less 

m 


than TrA 0 (X.) only if the f themselves decrease rapidly with increasing m. 
OCp m 


This, however, is not the case, for just as in Eq.(3.16), 
BZ 


Lim £ 

x.-<*> q 

l ~ 


uj(jq) i 


A, (q) 


X-f 


(A. 11) 


The origin of this behavior is the — behavior of • \ ■ as q- 

q (bUSA ~ 


«2°’ 43 


To circumvent this difficulty one must find a matrix M 0 (q) whose behavior 

ttp ~ 

1 

at the origin is the same as that of 2 — — — r- e (jq)e (iq). . , 

b i uj(jq) a £ {3 ~ * and which leads 

J /s/ 

. 3 

to an integral J d q ^ a g(q) cos(q*X) which can be evaluated analytically. 


BZ 


Then we write: 


fV ap (X l ) 


^BZh 
MN q 


Af 5 JT 35 T e a (ja) e e Uft) ) ■ M ap <fL) ] °' ,B( a‘Xd ) 


h BZ 

+ MN q M ap 


(q) cos (q* X . ) , 


(A. 12) 


and compute the first integral by the special point metDd. Since the integrand 


has no troublesome — behavior, its expansion coefficients f should then 

q m 

decrease rapidly, and the number of special points then needed for an accurate 

determination of A^„(X.) should be (and is in fact) correspondingly small. 

OCp 

To simplify the calculation, we have actually only treated the Trace 

of A (X.) in the above fashion, subtracting off a function M(q) whose 
ag ~ 

behavior as q -» 0 is approximately that of ^ . (As q -*0, 

~ o j UKjqj ~ 

£ ■ .. . . . — ■ ■ - x - ■ ■ , Where d(q) is a function of direction. \y e have approximated 

. UJ(jq) d(q)q ’ 
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d(q) with (S dfi — where the c .(q) are the three speeds of sound.) 
vj J q c (q) / J 

Tables A-I and A-II show the elements of A „ ( X ) , for X. in the first shells 

ag ~i 

(T = 0°K) , TrA n (X.), and A „ (0) for T = 0°K and T = 90°K. Three different 
CXp OCp 

(bcc) special point sets were used, with n = 8, 40, and 240. Although one 


can only expect TrA 


(X.) to converge well, the individual matrix elements 
y *°3- 


also show good convergence. 
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FIGURE CAPTIONS 


Fig. 1 Structure factor s(k) and the one-phonon contribution S_^(k) for Na at 

T=0°K along [100]. For all the structure factor and cross section plots, 
the ordinate is in absolute units. 

Fig. 2 Structure, factor S(k) for Na at T=0°K and T=90°K along [100]. 

Fig. 3 Structure factor .S(k) for Na at T=0°K and T=90°K along [110]. 

Fig. 4 Structure factor S (k) for Na at T=0°K and T=90°K along [110]. Note the 

expanded vertical scale. 

Fig. 5 Structure factor S(k) for Na at T=0°K and T=90°K along [111]. 

Fig. 6 Structure factor S(k) for Na at T=0°F: and T=90°K along [211]. 

Fig. 7 Structure factor S(k) for Na at T=0°K and T=90°K along [310]. 

6 7 

Fig. 8 Structure factor S (k) for Li and Li at T=0°K along [100]. 

6 7 

Fig. 9 Structure factor S(k) for Li and Li at T=0°K along [110], 

7 

Fig. 10 Structure factor S(k) for -.Li and Li at T=0°K along [111]. 

Fig. 11 Lines of equal value of the function l/c|q(k)J in a (001) plane of the 

2tt 

lattice reciprocal to the bcc lattice. R is the point — (0,0,0), P then 

3 . 

ott *5 i 2rr 

point — (— , — , 0), and S is the point — (3, 1, 0), where a is the lattice 
a 2 2 a 

constant. The numbers 1.00,0.50,0.33, and 0.25 indicate the relative value of the 
function. 

Fig. 12 Cross sections W(k) and W (k) for Na at T=0°K along [100]. 

Fig. 13 Cross sections W(k) and W (k) for Na at T=0°K along [110]. 

Fig. 14 Cross sections W(k) and W (k) for Na at T=0°K along [111]. 

Fig. 15 Cross sections W(k) and W (k) for Na at T=0°K along [100]. Note the ex- 

cl 

panded vertical and horizontal scales, and the position of k=2k^. 
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7 

Fig. 16 Cross section W(k) for Li at T=0°K along [110] for two different values 
of the core radius, 1^=1. 06 and r c =2.00. Note the .expanded horizontal 
scale. 

7 

Fig. 17 Cross section W(k) for Li at T=0 c K along [111] for two different values 
of the core radius, r c =1.06 and r c =2.00. Note the expanded vertical and 
horizontal scales. 

Fig. 18 Cross section with all Compton scattering subtracted, W(k)-S (k), for 

^Li at T=0°K along [111], for two different values of the core radius, 

r =1.06' and r =2.00. Note the expanded vertical and horizontal scales, 
c c r 



TABLE I. 


Structure factor of Sodium at T=0°K along [100]. S (k) vs. Q, where 
k = (1,0,0). 

El 



0 

S (k) 

0 

S CM 

0 

S(k) 

Q 

S (k) 

0,050 

3,34f>o4 

1,950 

4.71E-01 

3„6on 

?, 326-0 1 

5,250 

3.39E-01 

0,100 

6,75E*04 

1,970 

7,946-01 

3,650 

2 , 58£s0 i 

5,300 

3,516-01 

0,150 

1 ,03E*03 

1 ,9g0 

1 , 2 0 E 00 

3,700 

2,926-0} 

5,350 

3 .&56-0 1 

0,200 

1 , 4 l E-03 

1.998 

1.22E 0} 

3,750 

.3 « 4 1E*0 1 

5,400 

3,616-01 

0,250 

1 „biE-03 

2.00? 

1 , 22E 0} 

3,80 0 

4 , i 46-0} 

5*450 

3,996-01 

0,300 

2 ,256-03 

2,020 

1.24E 00 

3,850 

5 , 3 6 E * 0 1 

5,500 

4 , 20E-0 1 

0,350 

2,726-03 

2,030 

8 B 39£- 0 } 

3,900 

7,826-01 

5,550 

4g4ag-01 

0,000 

3,246-oS 

2,050 

S.16E-01 

3,950 

t,52E 00 

5,600 . 

4.756-0} 

0*450 

3 e SlE^OS 

2,100 

2.756-01 

3,970 

2 , 5 1 E Oo 

5,550 

5.14E-01 

0,500 ' 

4 e 42£-03 

2,150 

i ,976-01* 

3,980 

3.75E Oo" 

5,700' 

5,6bE-01 

0,550 

5,086-03 

2,200 

1,596-01 

3,998 

3 , 7 1 E 01 

5*750 

6,386-01 

0,600 

5 p 78K-03 

2,250 

1 , 376-° 1 

0,00? 

3.72E 01 

5.800 

7 • 45E-0 1 

0,650 

6 b 52E«q3 

2,3oo 

1.24F-01 

4,020 

3,806 OO 

5,850 

9,23E-01 

0,700 

7.31E-03 

2,350 

1 , lbE^O i 

4,030 

2,56£ 00 

5,400 

1,286 00 

0,750 

8 , 14E-03 

2,400 

t , 1 1E»0 1 

4,050 

1,57E 00 

5,950 

2,346 00 

Ot .600 

9,036-03 

2,450 

l,07E-0i 

“4,100 

g,40E«oi * 

5,970 

3,766 00 

0,850 

9,95E-03 

2,500 

l.OSE-Ol 

4,150 

6 , 0 0 £ « 0 i 

5,980 

5.526 00 

0,900 

U 10 E-J 02 

*2,550. 

1,04E“01 

4,200 

4 9 84E»0t 

5,998 

5 a 33E 01 

' 0,950 

l , 22E-02 

2,600 

1 B 04E-Ol 

4,250 

4, 17E-0J 

6,002 

5,336 01 

0,970 

i,26£-02 

2,650 

1.04E-01 

4,300 

3a 74E-0 1 

' 6 * o ?0 

5,54E 00 

0 o 9 8 0 

i«29E«02 

2,700 ' 

1.0UE-01 

4,350 

3 , 4bE*»0 1 

6 , 0 3 0 

3,776 00 

0,996 

l t 34£-02 

2,750 

1.05E-01 

4,400 

' 3.27E-0J 
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2.366 00 

1,002 

l, 356-02 

2 , 8 no 

1,066*01 
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3, 13E«0 1 
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1,316 00 
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1 ,406-02 
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1.08E-01 
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3,O4E*0! 
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9,696-01 

1,030 

1 ,43E-02 
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1,106-0} 
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2.97E-01 
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8,016-01 
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1,486-02 
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1, 13E-0I 
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2„93E«0l 
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7,046-01 

i. 100 ' 

t , b 4 6 — 0 2 
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1 ,146*01 
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2s 896-0 1 

6*300 

6,426-01 
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1 ,826-02 
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1,156-0} 
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2,876-0} 
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6,006-01 
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2,026-02 

2,998 

i, 166 - 0 } 
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2 , 86 E- 0 l 
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5,706-01 
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2 , 25£w02 
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1 , UE -01 
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2 , 868-01 
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5 , 4q£«Q 1 
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2,50E*»i)2 

3,0?0 

1 , l 8 E “0 1 
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2,676-0} 
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5,346-0} 
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2 b 79Ew02 

3,0 30 
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?g69e-oi 
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5.236-01 
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3, he-02 

3,050 

1,206-01 

4*950 

2„92E“0 J 

6,600 

S, 15E" n l 

1,450 

3,466-02 

3,100 

1,256-01 

4.970 

2 , 9 4 E » 0 1 

6,650 

5,096-01 

1,500 

3,9tE-02 

3,150 

1,316-01 

4 t 980 

2*95E-° 1 

6,700 

5 , 05E-0 1 

1,550 

4,436-02 

3,200 

1 ,376-01 

4,998 

2,976-01 

6,75° 

5, 02E-0 I 

1,600 

5. 076-02 

3,250 

1,446-01 

5,002 

2.97E-0 5 

6,800 

5 , 0 OE-0 } 

1,650 

5,896-02 

3,300 

1,526-01 

5,020 

2,996-0} 

6,850 

5,016-01 

1,700 

6 , 98E-02 

3,350 

1 ,816-01 

5,030 

3, OOE-O t 

6,900 

5 , 02 E -0 1 

1,750 

6,536-02 

3,400 

t . 7 IE-01 

5,050 

3,036-0} 

6,950 

5.05E-01 

1,800 

1 .09E-01 

3,4 5 o 

1,836-0} 

5,100 
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6,970 

5 , 07 E-O 1 

1,850 

1,486-01 

3,500 

1,96E”01 

5.150 

3,1 8E-0 } 
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5,086-01 

1,900 

2,286-0} 

3,550 

a^UE-oi 

5,200 

3,286-01 

6,998 

5,106*01 
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TABLE II. 


Structure factor of Sodium at T=0 °'k along [110]. S(k) vs. Q, where 
k= VQ (1,1,0). 

r*-j a. 


Q 

S ( k) 

Q 

S(k) 

0,025 

1 ,92E-04 

o,9?s 

2,8iE w 0l 

0,050 

3 a 8 6 g, 0 0 4 

0,985 

4 , 76f *> Q 1 

0,075 

5 p 83E»04 

0,990 

7.21E-01 

0,100 

7 ,84E-Qy 

0,999 

7.32E 00 

0,125 

9,9lE^04 

1,001 

7.35F. 00 

0,150 

1 » 2 1 E v 0 3 

1,010 

7,49E"0i 

0,175 

1 e 43E»o3 

1,015 

5, 04E«oi 

0,200 

1 ,67£=03 

1,025 

3,09E"01 

0,225 

l,92E-03 

1,050 

1 ,63E"01 

0,250 

2,1 6E-03 

1.075 

1.15E-0J 

0,275 

2,46E^03 

1,100 

9, UE«02 

0,300 

2,76E»03 

1 .125 

7,76E«p2 

0,325 

3,09E«o3 

1,150 

6 B 86F-o2 

0,350 

3 , 44t * 03 

1,175 

6 e 26E-02 

0,375 

3 ,62Ee>03 

t,2oo 

S,83E"02 

0.4QO 

4 ,23E*03 

1,225 

5,53E-02 

0,425 

4g68E«03 

1,250 

5,31E»02 

0,450 

5, 176^03 

1,275 

5, 17E*02 

0,475 

5 s 72E«03 

1,300 

5 „ 08E ra 02 

0,485 

5,95E«G3 

1,325 

5,03E«02 

0,490 

fe B 07E-03 

1 ,3.50 

5, 01E-02 

0,499 

6,29E®03 

1,375 

S.04E-02 

0,501 

6 , 34£«o3 

1,400 

5, 09E»02 
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6 ,57£«03 

1,425 

s, i7e«o£ 
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6,70E«03 

1,450 

5,20E®O2 

0 f 525 

6,96E»03 

1,475 

5.41E-02 

0,550 

7 , 72E=03 

1,485 

5.48E-02 

0,575 

8.56E-03 

1,490 

5,51E*’02 

0,600 

9,59E®03 

. 1,499 

5,57E«02 

0,625 

l ,06E«.02 

l , 5 0 i 

5,59E^02 

0,650 

1 , 1 8E-02 

1,510 

5,65E»0£ 

0,675 

1 ,33E-02 

1,515 

5,69E"02 

0,70° 

1 ,50E-02 

1,525 

5 B 78E™0g 

0,725 

U70E-02 

1,550 

6 „ 0 1 E"02 

0,750 

i,9SE«02 

1.575 

6,28E'-02 

0,775 

2 , 25Eo02 

1,600 

6, 60E-02 

0,80 0 

2 8 63E*02 

1,625 

6,98E«0? ‘ 

0,825 

3 „ 1 3 E « 0 2 

1.650 

7 t 4lE«02 

0,850 

3 e 6lE"0£ 

1,675 

7,93F.*02 

0,875 

4,76E»02 

1,700 

8 ,54E*»02 

0,900 

6, 19E-02 

1,725 

9,27E^02 

0,925 

8,6lE«02 

1,750 

1,026^01 

0,950 

1 ,35E*0l 

1,775 

1 , 1 3E**0 1 


0 

S (k) 

0 

S ( k) 

1,800 

1 „ 2 7 E " (> l 

2,625 

3 fi 74E"0l 

1,825 

t,45E«nt 

2,650 

1 ,828>01 

1 ,650 

1 , 69E=*0 1 

2,675 

1 .91E-01 

1,875 

2,03E-01 

2,700 

a e o3f>oi 

1,900 

2 , 5 4 E « 0 1 

2,725 

2, lfeF-Oi 

1,925 

3 , 39E-0 1 

2,750 

2.33E-Q1 

1,^50 

5 , 09E”0 l 

2,775 

a.soE-oi 

1 ,975 

1.02E 00 

2,800 

2.80E-01 

1,965 

1,70E 00 

2,825 

3, UE-01 

1 ,990 

2,56E 00 

2,850 

•3,590>Oi 

t ,999 

2 P 56E 01' 

2,875 

4 ,.2 IE "01 

2,001 

2.S6E- 0j 

2,900 

5, 14E-0J 

2,010 

2,60E 0Q 

2,925 

6,696-01 

2,015 

1.75E 00 

2,950 

9,7 6E**01 

2,025 

1.06E 00 

‘ 2,975 

t,90E 00 

2,050 

5,5Se°0i 

2,985 

3 , 1 2E 00 

2,075 

3,86fc>0i ' 

2,990 

4.65E 00 

2,100 

3.03E-01 

2,999 

4.60E 01 

2,125 

2 , 54{* n 0 1 ■ 

3 Tool 

4 , 60£ or 

2,150 

2 , 22E*»0 1 

- 3,oio 

4„69E 00 

2,175 

1 ,99E«0l 

3,015 

3 „ i 6E" 00 

2,200 

l , 83E«0 1 

3,025 

1 , 94 E 00 

2,225 

1 , 7 1 E>0 i 

3,050 

1,026 00 

2,250 

l,b3E-0l 

3,075 

7, 18E^0l 

2,275 

1 ,56E"0 1 

3,100 ' 

5 1 68E™'0 1 

2,300 

i .51E-01 

3,125 

4,79E~0l 

2,325 

1 ,48E»0i 

3,150 

4,20E"°1 

2,350 

1 .45E-01 

3,175 

3 , 79F- 0 1 

2,375 


3,200 

' 3,5oE“01 

2,4oo 

1 ,44g a Ql 

3,225 

3,27F“01 

2,425 

1,4 4E«Ci1 

3,250 

3,1 OF. -01 

2,450 

l ,aSE-0l 

3,275 

2.90F-O1 

2,475 

1 .47E-01 

3,300 

2 » 88E-0 1 

2,485 

1 1 40E“0 1 

3,325 

2.81E-01 

2,490 

1 ,48e»0l 

3,350 

2.76F-01 

2,499 

1 , 49E«0 1 

3,375 

2 , 7JE”0 1 

2,501 

t .49E-01 

3.400 

2.71E-01 

2,510 

1 , 5 0 E * 0 1 

3,425 

2 , 7OE-0 1 

2,515 

1 ,5 lE**0l 

3,«50 

2 e 7 1 E”0 1 

2,525 

1 , 5 2 E « 0 1 

3,475 

2 # 73F-01 

2,550 

1 , 5 6 E “ 0 t 

3.405 

2 , 7 4 E - 0 1 

2,575 

1 ,6lE*0l 

3,490 

2.74E-01 

2,600 

l ,67E“0t 

3,499 

2.75E-01 




IS 
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TABLE III. Structure factor of Sodium at T=0°IC along [111]. S (k) vs. Q, 
where k = — Q(l,l,l). 

cL 


» 

sc k ) 

0 

S(k) 

0 

SCk) 

fJ 

SCk) 

0,050 

4 S 51E«"0« 

1 ,950 

5 , 6 3 E " 0 } 

3,600 

4, 31 £-01 

5,250 

8,93E"01 

0,100 

9,316^04 

1,970 

9,248-01 

3,650 

4.36E-01 

5,300 

9 , 63E-0 1 

0,150 

1 ,4?e:-<o3 

1,960 

1.38E 00 

3,700 

4,528-01 

5.350 

9,628-0} 

0,200 

2,1 1E»03 

1,998 

1.36E 01 

3,750 

4 , 6 1 E «• 0 1 

5,400 

9,038-01 

0,250 

2*908-03 

2,002 

1.36E Oj 

3*800 

5,32E-0i 

5.450 

8,408-0} 

0,300 

3,86E-03 

2,0 20 

1,428 00 

3,650 

6,24E-oj 

5,500 

7.94E-01 

0,350 

5 , 1 4 1> 0 3 

2,030 

9.68F-01 

3,9oo 

6, 03E**0l 

5,550 

7,658-01 

0,400 

6*808-03 

2,050 

6 , 09E«0 1 

3,950 

1,30E 00 

5,600 

7.52E-01 

Q,4S0 

9* 1 IE-03 

2,100 

3,4Sf>01 

3,970 

1 , 95E 00 

5.650 

7.54E-01 

0,500 

1 * 24E-02 

2,150 

2 , 62E“0 1 

3,980 

2 * 75E 00 

5,700 

7 ,668-0 1 

0,550 

1.75E-02 

2,200 

2,268-01 

3.998 

2.44E 01 

5,750 

7.88E-01 

0,600 

2,50E»02 

2,250 

2, JOE-01 

4,002 

?,44e 01 

5,800 

8,288-01 

0,650 

3*47E>q2 

2,300 

2,058-01 

4,020 

2.76E 00 

5, P50 

8,96E-01 

0,700 

4*098-02 

2,350 

2,Q-9£-0 J 

4,030 

j ,96E 00 

5,900 

i,o?e 00 

0,750 

4 , 02E-02 

2,400 

2,208-01 

4,050 

1.32E 00 

5,950 

1,338 00 

0,800 

3, 73E-02 

2,450 

2.398-01 

4,100 

8,41 E-0 1 

5*970 

1,718 on 

0^850 

3®56E-02 

2,500 

2,7lE«0i 

4e 150 

6.79E-0J 

5,980 

2 , 18E 00 

0*900 

3,54£r02 

2,550 

3,198-01 • 

4,200 

6,04 E "01 

5,-996 

1,478 01 

0,950 

3,67£,02 

2,600 

3 ,85E-o 1 

4,250 

5,708-0} ' 

6,002 

1,478. 01 

0,970 

3*77M2 

2,650 

4,568-01 

4,300 

5 » 58E-0 1 

6,020 

2,168 on 

0,9«0 

3,83'E»02 

2,700 

4,76E«01 

4,350 

5,598-01 

6,030 

1.69E 00 

0,998 

3 e 96E-02 

2,750 

4 , 30E-Q i 

4,400 

5, 736"° 1 

6,050 

1,328 00 

1*002 

3 * 99E*»02 

2,800 

3 , 75£'"0 1 

4,450 

6 , 0 1 E « 0 J 

6,100 

1,038 00 

i p 0 a 0 

4,l5ep02 

2,650 

3 , 3 6 E " 0 \ 

4 , 5 0 0 

6 , 40E” 0 1 

6,150 

..9,208-01 

‘1,030 

4,25E«02 

2,90 0 

3 , 15E«0 i 

4,550 

7.15E-01 

6,200 

8,668 "01 

1,050 

4,47E>02 

2,950 

3,078-01 

4,600 

8,018-01 

6,250 

8,418-01 

1,100 

5,248-02 

2,970 

3,078-01 

4,650 

8*83E«01 

6,300 

8 « 32E-0 1 

1,150 

6e41&-02 

2,980 

3.088-01 

4 r 7 0 0 

8*978-01 

6,350 

8.328-01 

1,200 

8 t 2lE-02 

2,998 

3 , 1 1 E”0 1 

4,750 

8 , 32g" 0 1 

6,400 

8,418-01 

1,250 

1 *088-01 

3,002 

3,11F-01 

4,80 0 

7.55E-01 

6,450 

8,618-01 

1,300 

1 ,35£-o1 

3,02 0 

3 , tSE-o 1 

4,850 

6 , 998-0 1 

6*500 

‘ 8,938-01 

1,350 

1 ,42E-0t 

3,030 

3,1 £8-0 1 

4,900 

6e66E«01 

6.550 

9.35E-0J 

1 ,a00 

1 e 30Ew0l 

3*050 

3,25E«"01 

4,950 

6.54E-01 

6,600 

9,658-01 

1,450 

1 , 1 6E - 0 1 

3.100 

3.S1E-0! 

4,970 

6.S3E*0i 

6.650 

l.OgE 00 

1,500 

1 8 o7E»01 

3,150 

3 , 94E»sO 1 

4,980 

6(t53E"0 t 

6 , 7o 0 

1,038 00 

1,550 

1 .03E-01 

3,200 

4 S 60E”01 

4,998 

6,558-0}. 

6,750 

9,898-01 

1,600 

1 t oae-ol 

3,250 

5,485-01 

5,002 

6,568-05 

6,800 

9,458-01 

1,650 

1 ,098-01 

3,300 

6,298-0 j 

5,020 

6 , 5 9 E « 0 1 

6.850 

9,1 2E-01 

1 ,700 ' 

i.iee-oi 

3,350 

6 , 33E“0 1 

5,03 0 

6 , 62E*»0 1 

6,900 

8,918-03 

1,750 

1 ,33E-0l 

3,400 

5,71 F-01 

5,050 

6 e 69E“0 1 

6,950 

8,838-01 

1,600 

1 ,57E«01 

3,450 

5,06E-01 

5,100 

6..96E-01 

6,970 

8,838-01 

1 ,850 

2 # 0 1E-0 l 

3,500 

4 , 65E "0 1 

5,150 

7 , 42E"0 i 

6,980 

8,838-01 

1 ,900 

2,9 1 E«0 1 

3.550 

4,4oe«oi 

5,200 

B,10E«01 6,^9^ 

cgfo^ 

8,848-01 
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TABLE IV. 


Structure factor of Lithium-7 at T=0 °k along [111]. S(k) vs. Q, 
where k = — Q (1,1,1). 

3 . 


Q 

S(k) 

0 

S(k) 

Q 

S(k) 

a 

Stk) 

0 e 050 

1 # 02f>()3 

1,950 

1 , 021 ? 00 

3,600 

6,785*01 

5*250 

1,015 00 

0,100 

2 # 07E = 03 

1,970 

1.63E 00 

3,650 

6 , 86£<*0 1 

5,300 

1.04E 00 

0*150 

3,20E«»03 

1 ,9ft0 

2.39E 00 

3*700 

7,06E-01 

5,350 

1,045 00 

0*200 

4„46E«03 

1,998 

2 , She 01 

3,750 

7,395*01 

“5,400 

1,015 0 () 

0,250 

5,9iE«»o3 

2,002 

2 , 2 sf 01 

3,800 

7.93E-01 

5,450 

9,865 = 0 l 

0 , 30 ° 

7,68^03 

2, n 20 

2 , 43 E 00 

3,«30 

8,835=01 

5»500 

9,64E"0i 

0,550 

9 ,94E«>o3 

2,030 

1.68E 00 

3,900 

1 ,04£ 00 

5,550 

9,50E-Oi 

0,400 

1 t 30£>02 

2,050 . 

1 b 08£ on 

3,950 

1,405 00 

5,600 

9 , 445” 0 1 

0,450 

1,735=02 

2,100 

6,335=01 

3*970 

1,835 0O 

5,650 

9 t ase-oi 

0,500 

2,38E«02 

2,150 

4,64E«01 

3,960 

2.36E 00' 

5,700 

9,50E**0l 

0,550 

3 , 43E-02 

2 , 2 0 0 

4, UF.-01 

3,998 

1.65E 0j 

5,750 

9,595=01 

0,600 

5,t8E*02 

2*250 

3 , 82E"0 1 

4,002 

1.64E 01 

5,800 

9 „ 7 3E” 0 1 

0,650 

7„70E°02 

2,300 

3,708=01 

4,030 

2*35E 00 

5,850 

9,95E”01 

0 « ?Q0 

9 e 17E-02 

2,350 

3,725-01 

4,030 

1,605 00 

5,900 

1,035 00 

0,750 

8,44E>02 

2,400 

3,86E”0l 

4,050 

1,38F 00 

5,950 

1.08E 00 

0,800 

7 , 47E»0g 

2,450 

a, l6E«oi 

4,100 

1 , 05E 00 

5,970 

1,135 00 

0,850 

6,96E^02 

2,500 

4,66E-0i 

4,150 

9, 1 9 E = 0 J 

5,980 

1,205 00 

0,900 

6,85Ep02 

2,550 

5 ,44^ = 0 1 

4,200 

r.sse^o 1 

' 5,998 

2,915 00 

0,950 

7,09E«,02 

2*600 

6 . 56E"0 1 

4,250 

8,235=01 

fe , 0 0 2 

2,695 00 

0,970 

7*28E-02 

2', 650 

7,845-01 

4,300 

8, 1 g£*0 1 

6,020 

1 o 1 9E 00 

-0,960 

7.39E*02 

2,700 

8,145=01 

4,350 

8.13E-01 

6,030 

1,1 2E 00 

0,998 

7,6iE«-02 

2,750 

7, 16E-01 

4,400 

6,23E»01 

" '6,050 

1.07E 00 

1*002 

7„69E«Q2 

2,800 

6*245=01 

4,450 

8,465=01 

6,140 

1,025 on 

1,020 

7,99E-02 

2,850 

5*645-01 

4,500 

8 ,85E"0 1 

6,150 

9,975-01 

1,030 

8,1 8E-0 2 

2,900 

5 , 32E-Q 1 

4,550 

' 9*365-01 

6*200 

9,845-01 

1,050 

8 , 6 1 E » 0 2 

2.950 

5,235-01 

4,600 

9,955-01 

6,250 

9, 775=01 

1,100 

i t olE~0l 

2,970 

5 , 2 4 E ■* 0 1 

4,650 

1.05E 00 

6,300 

9,75E-01 

1,150 

1 p 24£«0t 

2,980 

5 * 265*=0 1 

4,700 

1,0SE 00 

6,350 

9.76E-01 

1,200 

1 ,6iE**0l 

2,908' 

5,29E«01 

4,750 

t , 0 0 e on 

6,400 

9,78E"0t 

1,250 

2,1 9 £>»o 1 

3,002 

5,30E*Oi 

4,600 

9 ,52E = 0 1 

6,450 

9,648-01 

1,300 

2q87E=01 

3,0?0 

5 , 35E-0 1 

4,850 

9 , 1 6 E 0 0 1 

6,500 

9,915-01 

1,350 

2 , 98E-Q 1 

3,030 

5,39E™01 

4,900 

8,965-01 

6,550 

1 , 00E 00 

1,400 

2 * 55E-0 1 

3,Or,o 

5,48E-0i 

4,950 

8,895-01 

6,600 

1*015 00 

1,450 

2,l8E“0l 

3,100 

5,825-01 

4,970 

r,89E«01 

6,650 

1,015 00 

1,50 0 

1 ,98E-01 

3,150 

6 , 38E«0 1 

4,980 

8.89E-0J 

6,700 

1,0 15 00 

1,550 

1 , 9yE»o 1 

3,200 

7,245-01 

4,998 

R ,9 1E®0 1 

6,750 

1 , 0 1 E 00 

1,600 

i,9ie»oi 

3,250 

8,4ie«01 

5*002 

8,915*01 

6,800 

9,995-01 

1*650 

2,02£«0 1 

3,3oo 

9,555-01 

5,020 

8,935-01 

6,850 

9.94F-01 

1,700 

2 , 22E«0 1 

3,350 

9 , 525-0 1 

5,030 

R.95E-01 

6,900 

9 , 9 0 E - 0 1 

1,75° 

2 , 5 2 E v 0 1 

3,400 

8*585*01 

5,050 

8,995-0 J 

6,950 

9.89E-01 

1,800 

3.02E-O1 

3,450 

7,755*0* 

5,100 

9.15E-01 

6,970 

9, 895-01 

1,850 

3,678=01 

3,500 

7, 195*01 

5,150 

9,405=01 

6,980 

9,695-01 

1,900 

5,53E*01 

3,550 

6,88E»01 

5,200 

9 . 735-01 

6,998 

9,895-01 
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2 o 

TABLE A-I. M° = %(2k f ) A° (in units of lO" ) . 

M(R) = %( 21 c 4 -) 2 • Tr A (R) (in units of 10 2 ). 
N is the number -of special points. 

(Na, T=0°K) . 



CO 

ii 

40 

240 

M d (T=0°K) 

3.4367 

3.4762 

3.4832 

M°(T=90°K) 

7.9897 

8.5890 

8.82 58 

M(g= (1,1,D) 

1.126 

1.134 

1.133 

M(R=(2,0,0)) 

0.538 

'0.541 

0.540 

M(R=(2j2 3 0)) 

0.283 

0.261 

0.259 

M(R=(3, 1, 1) ) 

0.240 

0.223 

0.221 

M(R=(2,2,2)) 

0.473 

0.479 

0.477 

M(R=(4,0,0)) 

0.174 

0.167 

0.164 

M(R=(3,3 3 1)) 

0.169 

0.152 

0.148 

M(R=(4,2,0)) 

0.140 

0.099 

0.095 

M(R=(4,2,2)) 

0.116 

0.137 

0.113 
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TABLE A-II . M_(R) = £(2kJ 2 A Q (R) (in units of 10~ 3 ) . 

ccp ~ ± ap ~ 

N is the number of special points. 

(Na, T=0°K) . 



N 

M 

M 

M 

M 

M 

M 



XX 

xy 

yy 

xz 

yz 

zz 

EL=( 1,1,1) 

8 

3.754 

2.610 

3.754 

2.610 

2.610 

3.754 


40 

3.780 

2.664 

3.780 

2.664 

2.664 

3.780 


240 

3.778 

2.666 

3.778 

2.666 

2.666 

3.778 

R=(2,0,0) 

8 

0.822 

0 

2.278 ' 

0 

• 0 

2.278 


40 

0.716 

0 

2.345 

0 

0 

2.345 • 

- 

240 

0.708 

0 

2.345 

0 

0 

2.345 

^(2,2,0) 

8 

1.278 

0.698 

1.278 

0 

0 

0.270 


40 

1.2-15 

0.740 

1.215 

0 

0 

0.184 


240 

1.207 

0:744 

1.207 

0 

0 

0.181 

R=(3, 1, 1) 

8 

0.715 

0.225 

0.842 

0.225 

0.316 

0.842 

r' 

40 

0.557 

0.230 

0.836 

0.230 

0.444 

0.836 


240 

0.541 

0.233 

0.832 

0.233 

0.448 

0.832 

R=(2,2,2) 

8 

1.578 

1.039 

1.578 

1.039 

1.039 

1.578 


40 

1.598 

1.133 

1.598 

1.133 

1.133 

1.598 


240 

1.589 

1.139 

1.589 

1.139 

1.139 

1.589 

^=(4,0,0) 

8 

0.581 

0 

0.581 

0 

0 

0.581 


40 

0.212 

0 

0.727 

0 

0 

0.727 


240 

0.186 

0 

0.730 

0 

0 

0.730 

R=(3,3,l) 

8 

0.680 

0.528 

0.680 

0.073 

0.073 

0.331 


40 

' 0.668 

0.458 

0.668 

0.120 

0.120 

0.188 


240 

0.653 

0.464 

0.653 

0.125 

0.125 

0.179 - 

3= (4, 2,0) 

8 

0.465 

0.275 

0.465 

0 

0 

0.465 

40 

0.356 

0.167 

0.400 

0 

0 

0.234 


240 

0.331 

0.171 

0.391 

0 

0 

0.227 

R=(4,2,2) 

8 

0.388 

0.304 

0.388 

0.304 

-0 

0.388 


40 

0.386 

0.208 

0.491 

0.208 

0.303 

0.491 


240 

• 0.363 

0.214 

0.482 

0.214 

0.314 

0.482 
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